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Abstract 


By using an asymptotic formula known for the numbers of Euler and Bernoulli it is possible to 
obtain an explicit expression of the nth digit of m in decimal or in binary, it also makes it possible to obtain 


the nth digit of 2”. The calculation is made from the two inequalities 
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By isolating 1 in both cases, we can derive an approximation of the latter. The approximation is so good 
that in fact it allows to extract the nth bit of 1, the nth decimal of m and even the nth decimal of dex", n € 
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Bernoulli numbers and 7t 


The 2 inequalities are taken from [3], p. 809. By isolating therefore, the number tin the 
first one obtains: 
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Here the B,, are positive and n is even. It is easy to verify that with n = 1000, the error on 
tm is of the order of 0.293193 x 107393, which is less than 271999. So the thousandth 
position of this expression is the 1000th bit of the number t. As soon as n = 10, the 
formula is valid: it allows to have the nth position. Note that given the parity of Bernoulli 
numbers, we obtain 2 bits at each n, that of rand n and that of rank n-1. This is valid from 
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10 because Bernoulli numbers are a little chaotic when n is small, the rapid growth of 
these numbers is felt from n = 12. 


The good question that comes next is: is it possible to calculate the nth decimal place? The 
answer is actually yes. Here's how. 


Given the known relationship with the Zeta function, it is possible to stretch the formula 
by adding enough terms to it to have an error of less than 10~”. In the book of Abramowitz 
and Stegun [3] precisely in the same chapter we have the formula 
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As we know, ¢(2n) can be represented as an infinite product due to Euler, by carefully 
taking the first 4 terms we therefore add them to our expression to have: 
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Knowing that the next term would be (1 — —) this assures us that the error committed 


will be less than 10’ (-n). Indeed, by checking with n = 1000 (and even 10000) we find an 
error of the order of 0.1271934403 x 1071°*3.We conclude that for each n in the 
expression it seems possible to obtain the nth decimal of the number Tt. The difficulty is 
then to be able to calculate B_2n when nis large. The current record for Bernoulli numbers 
is 100 million. So, we can calculate the 100,000,000th decimal of mt with this process. Note 
that the calculation of Bernoulli numbers can be done in several ways, one of which 
requires knowing tt with good precision. Until 2002, the Bernoulli calculation record used 
precisely this formula and some values of the series for ¢(n), see [5] on this subject. 
Currently, the calculation uses properties of congruences and the Chinese remainder 
theorem. 

The formula is precise enough to get even more if we raise to the power n. In fact we will 
have 


By checking with n = 1000 we obtain the first 1000 valid digits of 1*1000, so the 1000th 
position if we want. To get the 1000th position of x, just calculate 


position(n, x) = [10 {10”x }] 
2 


| | and { } being the integer and fractional parts. If calculating 1“n, the formula must be 
converted to treat the number as a string. 


Euler numbers and 1 


We take the expression of the book of A & S page 809, we isolate mt and we obtain: 
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Here the E,,, are the positive Euler numbers obtained by the series expansion of 1/cos(x). 
It is easy to verify that with n = 1000 we obtain a precision of 957 decimal digits. But 
knowing the error made and the first formula, we can deduce that this error is 
proportional to 9*(-1000). So, we can say that in base 8 the expression gives the nth digit 
of m in octal, a fortiori in binary. 

To go further as with the Bernoulli numbers, it is enough to add only 1 term to have an 
error lower than 10*(-n). Which gives us 
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And if we check with n >>1, for n = 1000 we have here an error on t of the order of 10’(- 
1198), so the nth decimal of m can be calculated with the formula using the position 
function (n,x). 

To get the nth digit in base 10 of z” using Euler numbers, we add terms until the error is 
smaller than 10™. We have to take more terms here since the powers of m take us away 
from the decimal point. So we finally have the following formula 
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For 1/ 1, in base 10 still with Euler, we can simply invert 
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It therefore makes it possible to obtain the nth decimal of 1/ 7 in base 10. The scope of 
this formula is limited by the ability to evaluate E_2n . With an Intel icore 9900K clocked 
at 5 ghz, I managed without too much difficulty to calculate E_(2,300,000). So position 
2,300,000 can be calculated by this process. 


In the end, we can therefore calculate the nth decimal digit of mand mn. In the end, we 
can therefore calculate the nth decimal digit of m and mn. The same process also makes 
it possible to calculate the nth decimal of the following numbers. 


e”™, In(z), Vx 
By isolating n! from the equation we obtain a much better formula than that of Stirling. 
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The same with Euler numbers 
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The 2 expressions are in absolute value and n even. A summary analysis of the error shows 
for n = 1000 shows that for the first the error is of the order of 271°, i.e. 10739 and 
371000 or 107-477 for the 2nd. With Stirling's formula, the error is of the order of 1072? only. 
As with the previous formulas, we can improve even more if we add terms drawn from 
the function ¢(n). 


But there is better, by comparing two successive terms of the expression with the 
Bernoulli numbers it allows to get rid of the exponent n. If we make the ratio, in fact we 
obtain for 77a closed expression. 


We had 
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By comparing 2 successive terms we get 
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Which is valid for base 4 given that the error is of the order of 4", by applying the same 
technique with Euler numbers we then have: 
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Whose error with respect to 1*2 is 9~”. So in base 9. We could use the same arguments 
for base 10 in both cases by adding terms. Hence the age-old question: is there an 


4 


interpretation of successive Bernoulli numbers allowing a combinatorial interpretation 
of the binary digits of 7, the question is plausible. For the first 20 terms, the thing is 
possible I believe. 


Real calculation of all the digits of 7 


The process would go further by calculating all the figures. Indeed, by pushing the formula 
with many more terms, we quite easily reach a few million digits, but the speed of 
convergence is not geometric, it is polynomial even if the degree is very high. If we take 
the largest known value of B_2n (100 million), we will therefore have the series (or 
product if we take the Euler product) as follows. 
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So, with n = 100 million, and even if we take 1000 billion terms, we will have at best 1.2 
billion digits of precision (1000 billion = 101%) so 12 digits only. This is a lot but little 
compared to the 62,000 billion digits obtained recently with the Chudnovsky formula, the 
latter has a convergence speed which is geometric, of the order of 14 decimal digits per 
term. 


Other asymptotic formulas containing the number tt 


First the factorial due to Stirling, 
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Which, once we isolate z and use the bootstrap method, 
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Will give some digits of z but even with n = 10? only 30 digits of precision is obtained. It 
is possible to stretch the first expression as desired but this convergence is polynomial, 
the gain is marginal. 


Another one used by Ramanujan and Hardy for the partition of n, p(n), is approximately 
equal to 
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Which also has a fairly slow convergence even if n = 10°, here we could extend the 
expression with Rademacher's formula, which is much more complex but which makes 


the number tT appear in several places, which makes the inversion operation for isolate 
rather difficult. This possibility cannot be ruled out at this time. 


With the central binomial coefficient, the following approximation holds. 
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Again, by isolating 7 we get 
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That can be pushed further by using the classic asymptotic expansion, 
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This time the approximation is the best so far, for n = 10° we get 100 digits of precision. 
The precision obtained if n >1 is much higher than the other 2 formulas but is not 
sufficient to obtain a precision of the order of the nth bit even if n is very large. If n = 10° 
we still obtain a hundred decimal places of precision. The asymptotic expansion is known 
and brings out the Bernoulli numbers which grow in size quite quickly beyond 14 terms. 
So the formulas obtained with Euler or Bernoulli numbers are by far the most accurate. 
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